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FOREWORD 

The  reseaich  described  in  this  repon  was  performed  during  the  1994  fiscal  year  as  part 
of  an  effort  to  improve  radar  classification  and  identification  capabilities  for  non- 
cooperative  airborne  targets.  This  problem  continues  to  be  a  primary  goal  of  radar 
research  programs  and  considerable  effort  has  been  expended  withm  the  last  few  decades 
in  attempts  to  solve  it.  The  accuracy  of  these  methods  depends  upon  available  radar 
resolution  and  so,  indirectly,  it  depends  on  radar  bandwidth. 

A  novel  technique  for  radar-based  non-cooperative  target  recognition  has  been  created 
and  verified  using  synthesized  data.  This  current  work  has  concentrated  on  verification  of 
the  method  using  measured  data  from  scale  model  targets  in  an  anechoic  chamber.  The 
method  uses  ordinary  azimuth  and  elevation  monopulse  tracking  data  and  is  appropriate 
for  application  in  many  existing  radar  systems. 
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INTRODUCTION 

Previously,  a  novel  technique  for  extending  the  information  available  from  range  profiles  of 
airborne  radar  targets  has  been  developed  and  applied  to  the  problem  of  automatic  target  recognition 
[References  1-4].  Since  so  much  of  the  target  shape  information  is  tied  up  in  the  croas>range  part  of 
the  radar  "image**,  the  principal  difficulty  with  down-range  profile  techniques  has  been  the  limited 
target  information  available  for  classification  purposes.  Usual  (e.g.,  Inverse  Synthetic  Aperture 
Radar,  or  ISAR)  schemes  are  often  unable  to  obtain  this  cross-range  information  because  accurate 
association  of  aperture  (induced  by  target  motion)  with  data  (collected  over  time)  is  very  difficult 
[see  Reference  3].  The  alternate  method  that  has  been  proposed  uses  data  collected  over  an  uaJcaown 
aperture  and  does  not  require  range-induced-phase  adjustments  across  the  set  of  measurements.  In 
addition,  the  technique  is  not  affected  by  aspect-induced  scintillation.  This  statistical  technique 
makes  use  of  the  fact  that  the  spatial  derivatives  of  the  scattered  field  phase  are  independent  of 
range-induced  phase  shifts  and  depend  upon  the  cross-range  extent  of  the  target. 

The  present  discussion  reports  on  initial  experimental  efforts  we  have  made  to  demonstrate 
the  feasibility  of  the  "enhanced  range  profile”  technique.  Using  scale  model  targets  in  an  anechoic 
chamber  we  have  obtained  very  encouraging  results  which  display  the  predicted  features  of  the 
algorithm  and  show  how  the  method  might  be  applied  to  automatic  target  recognition  (ATR)  and 
aimpoint  selection.  We  begin  by  reviewing  the  ^gorithm  itself  and  developing  the  relevant  ideas. 
Then  we  outline  the  measurement  process  and  show  some  of  our  results. 

BACKGROUND 

Phase  derivatives  have  been  used  for  target  tracking  purposes  for  many  years.  So-called  "phase 
monopulse”  systems  estimate  the  phase-gradient  of  the  scattered  field  by  differencing  the  outputs  of 
closely-spaced  antennas.  The  direction  of  this  gradient  will  lie  normal  to  surfaces  of  constant  phase 
(phase  fronts).  When  the  target  is  a  point  scatterer  the  phase  fronts  will  be  spherical  and  centered 
on  the  target,  and  their  normal  directions  will  indicate  target  bearing.  Of  course,  targets  are  uot 
individual  point  scatterers  and  the  phase-gradient  will  not  always  point  to  the  target.  The  difference 
between  actual  target  bearing  and  that  estimated  by  the  phase-gradient  is  known  as  target  "glint 
error”  and  depends  upon  overall  target  structure. 

A  short  radar  pulse  incident  on  a  target  can  be  used  to  isolate  separate  down-range  target 
elements  (see  Figure  1).  This  is  the  idea  behind  range  profiles,  which  are  maps  of  scatterer  reflectivity 
as  a  function  of  rwge.  By  combining  short  pulses  with  high  resolution  tracking,  target  subcomponent 
bearing  directions  can  sometimes  be  determined  as  a  function  of  range  as  well.  Unfortunately, 
targets  are  usually  sufficiently  complex  that  each  slice  will  generally  contain  more  than  one  point 
scatterer  and  scintillation  interference  effects  will  still  preclude  accurate  single-measurement  bearing 
estimates — rather,  statistical  estimates  based  on  multiple  measurements  are  required. 

When  a  set  of  phase  derivative  measurements  are  collected  over  a  short  time  interval  the  effective 
aperture  may  be  assumed  to  be  small.  If  the  target  aspect  cannot  be  accurately  determined  then 
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FIGURE  1.  Th«  RacUr  PuIm.  A  naiTo«  radar  pulie  incident  upon  a  target 
will  illuminate  only  a  range  ‘'•lice”  at  any  particular  inatant. 


the  glint  error  can  be  considered  to  be  a  function  of  random  aspect  angles.  Let  4  denote  the  phase 
of  the  scattered  field  so  that  Fteatt  =  |F«eatt|  exp  (t^).  Then  the  cross-range  component  of  the  phase 
gradient  along  a  direction  $  (see  Figure  2)  is  given  by 

ji(tf)  s  (icosfl +  58intf)  •  (1  ) 

This  will  be  the  ideal  measurement  made  by  two  closely  spaced  antennas  lying  along  the  0  direction. 

At  the  high  frequencies  employed  by  typicdi  phase  monopulse  “seekers”  the  phase  factor  4  will 
be  a  rapidly  varying  function  of  random  aspect  and  its  argument  can  be  modeled  as  a  uniform 
random  variable  on  (— ir,  t).  Let  p(x,  y)  denote  the  local  scattering  reflectivity  of  the  target  at 
position  («,]/)  (see  Figure  3).  When  the  target  is  complex  (the  number  of  scattering  centers  is 
large)  then  the  central  limit  theorem  implies  that  Fse»tt  »  normally  distributed  and,  because  of  the 
distribution  of  the  phase  factor,  has  zero  mean.  The  statistics  of  g  follow  from  a  straightforward, 
although  somewhat  lengthy,  calculation  [References  1,2,4].  The  probability  density  of  g  is  given  by 

“  2(«3  +  pJ)3/3’  ) 

where 

K*(fl)=  /  |p(«',y')|’(a:'co8^ +  p'sintf)*  dx'di/.  (3) 

Jwc> 

Equation  (3)  is  the  “second  electrical  moment”  of  the  slice  of  target  in  the  0  direction.  Target 
shape  information,  in  the  form  of  scattering  center  moments  weighted  by  loc^  scatterer  strength, 
can  be  deduced  from  a  set  of  measured  values  of  g  by  determining  the  parameter  that  best  fits 
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elevation  co(vd  axis 


FIGURE  3.  The  Rftder  CoordinAtee.  Hen,  f,Sf(fseO)  and  f«air(«s:v/3)  denote  the  components  of 
the  measund  target  bearing  in  the  asimuth  and  elevation  directions,  respectively. 
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the  atatiaticB  of  the  data  to  the  probability  density  of  Equation  (2).  Moreover,  it  is  easy  to  see  that 
K  (0)  =  \/k^  (6)  defines  an  ellipse  whose  shape  and  orientation  depend  upon  the  cross-range  extent  of 
the  target.  When  combined  with  time-domain  data  acquired  from  a  narrow  pulse,  Equation  (3)  will 
yield  a  family  of  down-range  ellipses  which  will  enhance  ordinary  down-range  profiles  by  providing 
additional  (although  limited)  cross-range  target  structure  information. 

The  data  used  to  estimate  the  “electrical  moments  of  inertia”  of  an  airborne  target  are  the 
ordinary  angle  of  arrival  (AOA)  measurements  used  in  radar  tracking  systems.  Since  an  airborne 
target  will  alter  its  orientation  over  time  as  it  maneuvers  (or,  even  if  it  is  flying  “strught  and  level”) 
the  AOA  will  fluctuate  with  time  as  well.  (In  the  straight-flight  case  the  target  orientation  will 
typically  vary  randomly  over  several  degrees  in  only  a  few  seconds  and  this  results  in  significant 
variation  in  the  AOA  estimates.)  The  statistics  of  this  variation  are  determined  by  the  target 
orientation  and  the  cross-range  target  structure  in  a  known  way.  By  fitting  the  statistics  of  the  AOA 
measurements,  collected  over  time,  to  this  model,  this  cross-range  information  can  be  determined 
for  each  slice. 

The  resulting  information  obtainable  is  overall  target  orientation  and  moment  “ellipse”  esti¬ 
mates  for  each  range  slice  (the  moment  ellipse  has  semi-axes  defined  by  the  electrical  moments  of 
inertia).  In  particular,  each  range  slice  is  mapped  to  an  ellipse  which  is  itself  defined  by  the  five 
parameters  described  by  Figure  4.  The  collection  of  all  range  slice  ellipses  associated  with  a  target 
at  a  particular  aspect  defines  the  total  available  target  structure  and  orientation  information. 


Elevation 


FIGURE  A.  EUipw  Parameter  Definition.  Theee  five  parameters  encode  the  croes-range 
target  information  available  for  each  down-range  eliee. 


Many  existing  radar  tracking  systems  are  capable  of  acquiring  this  data  with  only  minor  modi¬ 
fication.  For  our  purposes,  and  because  of  limited  funding,  we  chose  to  collect  sample  data  using  the 
system  described  in  Figure  5.  An  amplitude  monopulse  feed/comparator  was  attached  directly  to 
an  HP-8510  network  analyser  and  used  to  collect  frequency-domain  data  over  a  band  of  42-48  GHs 
from  a  scale  model  target  on  a  rotator  in  an  anechoic  chamber.  The  target  was  rotated  through 
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2  degrees  in  azimuth  and  3  degrees  in  elevation  and  was  numerically  transformed  (DFT)  into  the 
time-domain  as  part  of  the  preprocessing.  Finally,  50  of  these  time-domain  AO  A  measurement  sets 
(corresponding  to  50  randomly  selected  aspects  from  within  the  2  degrees  x  3  degrees  measurement 
domain)  were  fed  into  the  parameter  fitting  algorithm — the  results  from  the  target  of  Figure  6  are 
plotted  in  Figures  7-9. 


FIGURE  5.  Angle  of  Arrival  Measurement  System. 
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FIGURE  7.  Enhanced  Down-Range  Profiles  of  a  Scale  Model  Target.  Shown  are  the  cross- 
range  semi-axes  a  and  t  of  the  “inertia”  ellipses  plotted  as  a  fimctiou  of  range. 


FIGURE  S.  The  Ellipse  Center  Coordinates  r,  and  Plotted  as  a  Faction  of  Range. 
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FIGURE  9.  Th«>  EUipM  Orientation  Angle  a. 


DISCUSSION 


Because  the  anechoic  chamber  was  small  aod  precluded  the  use  of  large  scale  model  targets, 
the  conclusions  that  can  be  drawn  from  these  results  are  limited.  The  target  represented  in  Figure 
6  was  a  1/56  scale  model  Tu-22  (“Blinder”)  bomber  and  the  radar  system’s  center  frequency  and 
bandwidth  scale  to  800  MHz  and  ICO  MHz,  respectively.  Consequently,  the  results  are  not  nearly 
as  good  as  should  be  expected  from  most  currently  deployed  radar  systems  and  full  size  targets. 

To  better  appreciate  these  results  we  have  included  a  high  resolution  ISAR  image  of  the  same 
scale  model  target  using  comparable  data  (in  frequency  and  bandwidth)  (Figure  6).  This  image 
defines  the  electrical  target  and  contains  all  of  the  information  available  to  the  radar.  Note  that  the 
target  is  electrically  represented  by  less  than  about  20  scattering  centers  at  the  frequencies  employed. 
This  severely  taxes  the  assumptions  of  the  algorithm.  In  particular; 

—  The  estimate  of  ««  u  not  constant  along  the  target  axis  (as  it  should  be  idedly).  This  is  because 
there  were  too  few  scattering  centers  to  achieve  a  “good”  average.  Moreover,  these  errors  also 
^ect  the  semi-major  axis  estimates.  Higher  effective  frequencies  (as  are  appropriate  to  actual 
radar  systems)  will  result  in  more  effective  scattering  centers  and  so  increased  accuracy  in 
parameter  estimates. 

—  The  target  was  dominated  by  the  engine  inlets.  Since  the  cross-range  ellipses  are  determined 
by  intensity  weighted  averaga  of  the  scattering  center  position,  the  ellipse  associated  with  this 
portion  of  the  target  is  much  more  narrow  that  would  be  expected  by  a  casual  inspection  of  the 
target  model. 
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In  fact,  the  overall  target  is  electrically  much  more  narrow  than  it  is  physically.  (Note  that  the  wings 
an  not  electrically  active  beyond  the  pods.)  Since  so  much  of  the  strong  features  are  concentrated 
along  the  center  of  the  target,  this  makes  it  difRcult  to  interpret  the  data  results,  which  contain  such 
small  variation. 

Nevertheless,  the  results  of  Figures  7-9  clearly  demonstrate  the  effectiveness  of  the  algorithm 
even  for  this  worst-case  test.  The  target  orientation  has  been  accurately  estimated  as  well  as  the 
overall  target  shape  and  structure. 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  do  not  feel  comfortable  claiming  the  algorithm  as  “proved  effective  for  ATR."*  Rather, 
because  of  the  limitations  in  the  test  we  have  only  demonstrated  that  the  algorithm  is  remarkably 
effective  even  when  the  model  assumptions  are  significantly  relaxed.  In  particular,  the  technique 
is  clearly  appropriate  for  the  more  limited  problem  of  aim-point  selection  in  rf-guided  missiles.  In 
fact,  because  the  identical  data  are  also  used  for  guidance,  this  method  may  prove  to  be  an  optimal 
(practical)  aim-point  selection  algorithm. 

A  more  complete  demonstration  of  the  method  (for  ATR)  will  require  larger  model  targets 
with  more  realistic  complex  structures.  However,  given  the  relative  success  so>far,  as  well  as  the 
experience  gained  in  the  current  effort,  a  complete  and  convincing  demonstration  should  be  straight 
forward. 
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